The Clustering Dipole of the Local Universe from the Two 
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ABSTRACT 

The unprecedented sky coverage and photometric uniformity of the Two Micron All Sky Sur- 
vey (2MASS) provides a rich resource for investigating the galaxies populating the local Universe. 
A full characterization of the large-scale clustering distribution is important for theoretical studies 
of structure formation. 2MASS offers an all-sky view of the local galaxy population at 2.15/um, 
unbiased by young stellar light and minimally affected by dust. We use 2MASS to map the 
local distribution of galaxies, identifying the largest structures in the nearby universe. The inho- 
mogeneity of these structures causes an acceleration on the Local Group of galaxies, which can 
be seen in the dipole of the Cosmic Microwave Background (CMB). We find that the direction 
of the 2MASS clustering dipole is 11° from the CMB dipole, confirming that the local galaxy 
distribution accelerates the Local Group. From the magnitude of the dipole we find a value of 
the linear bias parameter b = 1.37 ± 0.3 in the K s -ba,nd. The 2MASS clustering dipole is 19° 
from the latest measurement of the dipole using galaxies detected by the Infrared Astronomical 
Satellite (IRAS) suggesting that bias may be non-linear in some wavebands. 



Subject headings: galaxies:clusters 

1. Introduction 

Since the discovery of a dipole in the CMB (Lu- 
bin et al. 1983; Fixsen et al. 1983) it has been 
identified as the Doppler signature of our motion 
relative to the CMB rest frame. Strong evidence 
in favor of this interpretation has been provided 
by the observation of a similar dipole in the sur- 
face brightness of radio galaxies in the same di- 
rection (Blake & Wall 2002). Blake & Wall make 
the distinction between velocity dipoles caused by 
our motion and clustering dipoles which measure 
the distribution of galaxies. Attempts to mea- 
sure the clustering dipole of local galaxies in the 
1980s (Meiksin & Davis 1986; Yahil et al. 1986; 
Lahav 1987; Villumsen & Strauss 1987) from only 
the fluxes and positions of galaxies led to re- 
sults that were within 30° from the CMB velocity 
dipole. The inclusion of redshift information led 
to a controversy over out to which distance the 
clustering dipole converges (Strauss & Davis 1988; 
Lynden-Bell et al. 1989; Rowan-Robinson et al. 



1990; Plionis & Valdarnini 1991; Scaramella et al. 

1991) . Through the 1990s progressively larger red- 
shift surveys were used in an effort to determine 
where the clustering dipole converges and to im- 
prove the calculation of cosmological parameters 
(Strauss et al. 1992; Hudson 1993; Schmoldt et al. 
1999; Rowan-Robinson et al. 2000). These im- 
provements included determining the shot noise 
in the sample, optimizing the window function for 
smoothing and estimating the contribution of non- 
linear effects to the velocity dipole (Strauss et al. 

1992) . These efforts have recently determined that 
the clustering dipole converges at 150 — 200/i _1 
Mpc (Schmoldt et al. 1999; Rowan-Robinson et al. 
2000). 

2MASS (Skrutskie et al. 1997) is the first 
near-infrared (JHK S passbands) all-sky survey. 
2MASS has an effective image resolution of 1 arc- 
second and 100 x greater sensitivity than the far- 
infrared all-sky survey, IRAS. Most passbands 
tend to be sensitive to star formation rate, while 
the K s passband is most sensitive to stellar mass 
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(Bell & de Jong 2000), which probably makes the 
-RT s -band a better tracer of total mass. Therefore, 
2MASS offers a unique opportunity for calculating 
a flux-weighted clustering dipolc. 

The median depth of the survey is z — 0.073, 
or 220ft,- 1 Mpc (Bell et al. 2003), a distance past 
where the dipole has been shown to converge. 
Figure f shows the distance distribution for a 
magnitude-limited sample of 2MASS galaxies as 
determined from 6, 104 redshifts from the Sloan 
Digital Sky Survey early data release (Stoughton 
et al. 2002). At a limiting magnitude of K s = 
13.57 mag, 2MASS includes all and brighter 
galaxies out to 200h~ 1 Mpc. 

We can use the 2MASS luminosity function 
measured using the Sloan redshifts (Bell et al. 
2003) to determine the completeness of the sur- 
vey as a function of distance. The result, shown 
in Figure 1, is that the survey contains 60% of the 
total luminosity and 76% of the total flux inside 
200/i _1 Mpc. Thus 2MASS is a fairly complete 
survey of the total flux in the local universe. As 
we expect the clustering dipole to be dominated 
by bright objects, any missing flux is unlikely to 
contribute significantly to the dipole. Later we 
show that this is true by demonstrating the con- 
vergence of the clustering dipole as a function of 
apparent magnitude. 

2. The 2MASS Catalog 

The 2MASS final release extended source cata- 
log is the most complete resource of local galaxy 
observations at near-infrared wavelengths. The 
catalog contains over 1.6 million sources, and for 
most Galactic latitudes (|b| > 20°) has > 90% 
completeness and > 98% reliability down to K a < 
13.5 mag (Jarrett et al. 2000b). For 5 < |6| < 20°, 
in regions where the stellar density is less than 
10 4 per square degree with K s < 14 mag stars, the 
catalog maintains its high completeness and is still 
> 80% reliable (Jarrett et al. 2000a). Hereafter, 
we adopt the Kron (1980) magnitudes, which at- 
tempt to recover the "total" galaxy flux using 
apertures related to the galaxy radius (limited in 
2MASS to a 5"minimum). 2MASS Kron magni- 
tudes underestimate the true total flux systemat- 
ically by 0.1 mag for K s > 11 (Bell et al. 2003), 
which we correct for when converting magnitudes 
to flux. The extended source catalog is 97.5% 
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Fig. 1. — Properties of the 2MASS galaxies as 
a function of distance R. The top panel shows 
the percentage of 2MASS galaxies in 20ft~ 1 Mpc 
bins as a function of distance. The middle panel 
shows the percentage of luminosity that falls be- 
low our flux limit of K a = 13.57 mag as a function 
of distance. The bottom panel shows the percent- 
age of the total integrated flux out to a given dis- 
tance that is above the flux limit of K s — 13.57 
mag (solid line). Also shown is the percentage of 
the total integrated luminosity (dashed line) as a 
function of distance, where both calculations have 
used the measured luminosity function (Bell et al. 
2003). 

complete within the Sloan early data release for 
extinction-corrected Kron K s < 13.57 mag (Bell 
et al. 2003). While the reliability decreases at 
lower galactic latitudes we will adopt this limiting 
flux and explore the effect of a brighter cutoff later 
on. The photometry of the catalog is very uniform 
with errors of less then 0.1 magnitudes (Nikolaev 
et al. 2000). We have performed an analysis of the 
galaxy density with a number of possible contami- 
nants and find that the cross-correlation is always 
less then 5% (Mailer et al. 2003). 

We apply an extinction correction based on 
Galactic dust maps (Schlegel et al. 1998) to cre- 
ate a foreground corrected catalog of galaxies with 
K s < 13.57 mag. From the overlap with the 
6, 104 Sloan spectroscopically confirmed galaxies 
we find that bright (K s < 12 mag) galaxies have 



2 



dust-corrected colors narrowly clustered around 
0.8 < J — K < 1.2. Thus, we exclude objects 
brighter than K s = 12 mag with J — K < 0.6 or 
J — K > 1.4, as there are no galaxies with these 
colors found in the confirmed sample (~ 2000 ob- 
jects). Since the photometric errors increase the 
spread in J — K colors for faint objects we do not 
extend this color criteria below K s = 12. 

In the region I > 230° and I < 130° we mask out 
the region \b\ < 7° and extend this to |fr| < 12° in 
the region I > 330° and I < 30°, where the 2MASS 
sources are clearly obscured. We also remove 17 
galaxies that are members of the Local Group. 
This leaves us with 748, 080 galaxies covering 90% 
of the sky. Note that applying no dust correction 
reduces the number of galaxies to 704, 035. 

We display the most comprehensive represen- 
tation of the local Universe in Figure 2. We see 
a complex web of galaxies stretched out along fil- 
aments between dense nodes often corresponding 
to the location of giant superclusters. The most 
obvious structure is the "Great Attractor" found 
within the region bounded by 300 < I < 360° 
and —45 < b < +45°. This structure is roughly 
divided into north and south Galactic hemisphere 
complexes - the overlapping Hydra-Centaurus and 
"Shapley" supercluster complex (0 < b < +45°) 
and the Pavo-Indus wall (—45 < b < 0°) (Dressier 
1988; Fisher ct al. 1995). The Perseus-Pisces 
"chain" extends from the Perseus supercluster (I = 
150°,6 = -15°) to (I = 110°,& = -35°)(Haynes 
& Giovanelli 1986), and appears to reach up to 
b = +30° along latitude 160° (Fisher et al. 1995). 
It has been demonstrated that there is no phys- 
ical link between the Perseus-Pisces chain and 
the Pavo-Indus wall(di Nella et al. 1996). Addi- 
tional features include the Horologium-Reticulum 
supercluster (/ = 265°, b = —55°) (Lucey et al. 
1983), Virgo supercluster (I = 285°,6 = +75°), 
the Ursa Major cloud (I = 145°,6 = +65°), 
and the NGC 1600 group (I = 210°,6 = -30°) 
(Fisher et al. 1995). Finally, the "Great Wall" at 
cz ~ 8000 km/s (Geller & Huchra 1989) is the 
faint feature running from the Hercules superclus- 
ter (I = 30°,6 = +45°) to (I = 180°,6 = +30°). 

3. Measuring the Clustering Dipole 

The inhomogeneity of these structures causes a 
gravitational acceleration on the Local Group of 



galaxies. The CMB velocity dipole (Lincweaver 
et al. 1996) results from the motion of the Sun 
relative to the CMB standard of rest. This motion 
can be broken down into two parts; the motion 
of the sun relative to the Local Group and the 
motion of the Local Group with respect to the 
CMB. Using the most recent values (Courteau & 
van den Bergh 1999) we find that the Local Group 
velocity relative to the CMB is 622 kms -1 in the 
direction l cm b = 272, b cm b — 28. 

Since both gravitational force and flux fall off as 
distance squared, the net acceleration of the Local 
Group is proportional to the dipole of the light 
distribution for a constant mass-to- light ratio T. 
Thus, 



(1) 



where Si is the flux received from each galaxy and 
fj is a directional unit vector. However, the light 
may be a biased tracer of mass (Kaiser 1984); if 
there is linear biasing then the true acceleration 
would be a factor 1/6 times the measured value. 

Linear theory can be used to estimate the resul- 
tant velocity (Peebles 1980), which depends on the 
matter density of the universe Cl m . The expected 
velocity of the Local Group is therefore, 



v = 



2 f(Slr. 



3 bH n r . 



(2) 



where f{fl m ) ~ fl° m 6 (Peebles 1980) and H is the 
Hubble constant. Of course, this is only true in 
linear theory. Nonlinear N-body simulations show 
that the velocity of a region like the local group is 
typically within 7° of the computed linear acceler- 
ation and that the magnitude agrees to within 20% 
(Davis et al. 1991). Measuring the light dipole 
therefore achieves two aims: (1) It verifies that 
the CMB dipole is truly caused by the Sun's mo- 
tion and (2) it can determine a combination of T^, 
Ct m and b. Since Q m and Tk can be measured by 
other means we can therefore determine the value 
of b. 

To accurately calculate a clustering dipole we 
need to consider the effect of the masked region. 
We address this in two ways: by cloning the sky 
above and below the masked regions (Lynden-Bcll 
et al. 1989), and by filling the masked region with 
randomly chosen galaxies such that it has the same 
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Fig. 2. — The upper two images show the southern and northern Galactic hemispheres in Lambert's projec- 
tion. The bottom image shows the whole sky in an Aitoff projection. The color denotes the smoothed flux 
density measured in units of the solar _?G-band flux per square degree. In the lower image the Galactic center 
has been masked out (dashed line) and adjacent regions have been cloned to fill the masked region. We have 
labelled the most prominent structures: (S) "Shapley" complex, which includes Hydra-Centaurus in the 
foreground; (P-I) Pavo-Indus wall; (P-P) Perseus-Pisces chain; (H-R) Horologium-Reticulum supercluster; 
(V) Virgo supercluster; (UM) Ursa Major cloud; (N) NGC 1600 group; and (GW) the "Great Wall" . The 
green line delineates the super galactic plane along which many of these structures lie. Also the CMB (white 
cross) and 2MASS (green cross) dipoles are shown. 
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Fig. 3. — The growth of the clustering dipole is 
shown as a function of the number of galaxies 
used to compute it ordered by their flux. The 
corresponding limiting K s magnitude is shown on 
the top axis. The dipole quickly rises to 80% of 
its final value. The faintest 300, 000 galaxies only 
change its value by < 5%. The growth of the total 
flux (dashed line) for the same galaxies is shown 
for comparison. The total flux rises more slowly 
and continues to rise even for the faintest magni- 
tudes in our sample. 



surface density as the unmasked area. We ap- 
ply both methods to estimate the uncertainty con- 
tributed by the masked region. 

Figure 3 shows the convergence of the magni- 
tude of the clustering dipole as a function of the 
flux limit of the survey. The dipole mostly arises 
from galaxies with K s < 12 mag, the brightest 
50, 000 galaxies. The dipole converges at fainter 
magnitudes and the addition of 300,000 galaxies 
with K s > 13.2 changes the dipole by less then 
5%. For comparison, the dashed line shows the 
growth in the total flux as a function of the num- 
ber of galaxies used, which continues to grow for 
faint magnitudes. This is compelling evidence that 
the dipole has converged in the 2MASS sample 
and that any flux we are missing in low luminos- 
ity galaxies is a minor contribution. 

To estimate the effect of shot noise on our cal- 
culation we perform bootstrap resampling on the 
galaxy catalog. We first bin the galaxies by their 
fluxes and then resample each flux bin to ensure 
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Fig. 4— The location of the 2MASS cluster- 
ing dipole is shown relative to the CMB velocity 
dipole. Also shown is the location of the cluster- 
ing dipole calculated from the IRAS Pscz catalog 
(Rowan- Robinson et al. 2000) . The position of the 
dipole with the masked area filled by cloning (tri- 
angle), filled randomly (diamond), or including no 
dust correction (square) are also shown. The error 
bars reflect the systematic uncertainties caused by 
dust and the mask. 



that the resampled catalog has the correct flux 
distribution. Performing the resampling 100 times 
the standard deviations of the clustering dipole di- 
rection is = 0.5°, Ab = 0.3° and in magnitude 
it is 0.5%. We find that the systematic uncertain- 
ties are much larger than the shot noise. 

The most important sources of systematic er- 
ror in our calculation are our treatment of the 
mask and of dust corrections. Cloning the sky 
above and below the mask gives a dipole point- 
ing towards I = 280°, b = 39°. If instead of 
cloning the adjacent sky we fill the mask with ran- 
domly selected galaxies the direction of the dipole 
changes to I = 277°, b = 37°. If we perform no 
dust correction on the galaxy catalog then the 
dipole moves to I = 277°, b = 40° (Figure 4). We 
adopt the centroid of these three measurements 
kipoie = 278 ±2.5°, b dipo i e = 38 ±2° as the best fit 
dipole with error bars reflecting the dust correc- 
tion and mask filling uncertainties. The dipole cal- 
culated only using galaxies brighter then K s = 13 
points to / = 278°, b = 40° verifying that the miss- 
ing fainter galaxies are not introducing any signif- 
icant error. 

The clustering dipole is 11° from the CMB ve- 
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locity dipolc, closer then any previous determina- 
tion. This is somewhat larger then the average 
separation of 7° found in the N-body simulations, 
but within 95% confidence limits. Therefore we 
conclude that the 2MASS determined clustering 
dipole is consistent with the direction of the Local 
Group motion given the uncertainties caused by 
non-linear effects and systematic errors. 

Interestingly our clustering dipole is 19° from 
the dipole measured most recently using the IRAS 
Pscz catalog (Rowan- Robinson et al. 2000). This 
is larger than the systematic errors in cither cal- 
culation and therefore suggests differences either 
caused by the passband of the survey or in the 
methods used to calculate the clustering dipole. 
It is possible that biasing schemes more compli- 
cated then simple linear biasing may be needed to 
explain this disagreement. The magnitude of the 
2MASS clustering dipole is 1.25 ± 0.06 x 1O" 3 L 
kpc~ 2 , where again the errors are dominated by 
the systematic uncertainties in the dust correction 
and the treatment of the masked area. 

4. Estimating b 

The Wilkinson Microwave Anisotropy Probe 
(Bennett et al. 2003) has recently determined the 
matter density of the universe to be TO = 0.27 ± 
0.04 and H = 71 ± 4kms~ 1 Mpc~ 1 (Spergel et al. 
2003). From the 2MASS luminosity function we 
can determine the K s luminosity density of the 
universe to be j K = 5.27 ±0.11 x 10 8 hL Q Mpc^ 3 
(Bell et al. 2003) and therefore the average mass- 
to-light ratio, Tx = 71ft. -1 ± 5L©/M Q . Using 
these values in eq. 2 we find that the expected ve- 
locity in linear theory is 8566 -1 ± 80 kms" 1 com- 
pared to the value measured from the CMB of 
622 kms -1 . Including the 20% uncertainty be- 
tween the linear prediction and the measured ve- 
locity in nonlinear N-body simulations we find 
b = 1.37 ± .34. This can be compared with a 
measurement of bias from the 2df Galaxy Red- 
shift Survey using rcdshift space distortions which 
finds that b = .94 ± .15 in the 6j-band (Peacock 
et al. 2001). 

If bias is luminosity or color dependent then 
we would expect a higher value from 2MASS 
than from the 2df Galaxy Redshift Survey. Since 
the measurement of bias function of wave- 
band and scale is one of the most important tools 



used to understand the formation and evolution of 
galaxies, this measurement of the clustering dipole 
not only confirms basic cosmological predictions, 
but will also lead to a better understanding of how 
galaxies form. 
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